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Abstract 

At very high energies scattering amplitudes in a spontaneously broken gauge theory into 
multi-particle final states are known to grow factorially with the number of particles pro¬ 
duced. Using simple scalar field theory models with and without the VEV, we compute total 
cross-sections with up to 7 particles in the final state at the leading order in perturbation the¬ 
ory with MadGraph. By exploring the known scaling properties of the multi-particle rates 
with the number of particles, we determine from these the general n-point cross-sections in 
the large-n limit. In the high-multiplicity regime we are considering, n^> 1 and An=fixed, 
the perturbation theory becomes strongly coupled with the higher-order loop effects con¬ 
tributing increasing powers of A n. In the approximation where only the leading loop effects 
are included, we show that the corresponding perturbative cross-sections grow exponentially 
and ultimately violate perturbative unitarity. This occurs at surprisingly low energy scales 
~ 40 — 50 TeV with multiplicities above ~ 150. It is expected that a repair mechanism or 
an extension of the theory has to set-in before these scales are reached, possibly involving a 
novel non-perturbative dynamics in the a priori weakly coupled theory. 



1 Introduction 


We are interested in scattering processes at very high energies into n-particle final states in the 
limit ?i> 1. In this case the well-known problem of divergences of large orders of perturbation 
theory [1, 2, 3, 4], is realised instead at the leading order. This is because even the leading- 
order Born diagrams for the n-point scattering amplitudes are expressed in terms of Feynman 
diagrams with large numbers of vertices, and with the numbers of diagrams growing factorially 
with n. At sufficiently high energies the production of high multiplicity final states, with n 
greater than the inverse coupling constant, is kinematically allowed and the n-point scattering 
amplitudes near the multi-particle mass thresholds grow as n! at leading order - i.e. tree level 
in a weakly coupled theory. 

In the simplest scenarios, production rates for such final states can be considered in a 
quantum field theory with a single scalar field of mass M and the coupling A. The model with 
a non-vanishing VEV ( h) = v, 

£ ssb = l d r hdiih _ ±(h 2 -v 2 ) 2 , (1.1) 

is a simplified version of the Higgs sector of the SM in the unitary gauge, describing neutral 
Higgs bosons of mass M = y/2Xv. We will refer to this model as the theory with spontaneous 
symmetry breaking (SSB) and it will be our principal case of interest for this paper. In addition 
to (1.1) will also consider multi-particle amplitudes in an even simpler theory, with no 
spontaneous symmetry breaking, 

£ noSSB = 1 W) 2 _ l _ M 2^2 _ ^4 (1 .2) 

The model (1.2) with an unbroken Z 2 symmetry has been widely used in computations of 
multi-particle rates in the 1990’s as reviewed in Ref. [5] and other papers referred therein and 
below. 

Our goal here is to compute the multi-particle rates directly in perturbation theory using 
one of the current state of the art publicly available numerical techniques, in the current case 
- MadGraph 5 [6]. The continuation procedure from moderate values of n = 7 particles in 
the final state, where our calculations are performed, to the regime with n 10 2 - 10 3 will be 
set-up and carried out in Section 3 based on the known scaling properties of the multi-particle 
cross-sections with n, as will be outlined next in Section 2 . 


2 Multi-particle production rates and the holy grail function 

Let us consider the multi-particle limit n 2> 1 for the n-particle final states, and scale the 
energy y/s = E linearly with n, E oc n, keeping the coupling constant small at the same time, 
A oc 1/n. It was first argued in [7] (for a review of subsequent developments see [5]) that in this 
double-scaling limit the production cross-sections a n have a characteristic exponential form, 

a n ~ e nF ( An > £ ) ) for n —>• 00 , An = fixed, e = fixed , (2.1) 


1 



where e is the average kinetic energy per particle per mass in the final state, 


e = (E — nM)/(nM), (2.2) 

and F(Xn, e) is a certain a priori unknown function of two arguments. F(Xn, e) is often referred 
to as the ‘holy grail’ function for the multi-particle production 1 in the perturbative sector. 

At small values of e and An, the large-n behaviour in (2.1) has been verified explicitly and the 
function F(Xn,e ) was computed in [7] for the VEV-less scalar theory (1.2), and later in [8] for 
the theory with the VEV (1.2) and more generally in a Gauge-Higgs theory. These computations 
were carried out in perturbation theory at tree-level combined with the simplifications arising 
in the non-relativistic limit e <C 1 for the final state particles. This approach has allowed for the 
analytic derivation the corresponding tree-level amplitudes and their phase-space integration 
for all values of n > 1. It was found that the dependence of the holy grail function on its two 
arguments, A n and e, factorises into individual functions of each argument 

logcr* ree In-Kxj nF tree (Xn, e) = n(f 0 (Xn) + /(e)) , (2.3) 


and the two independent functions are given by the following expressions in the model (1.2) 
Ref. [7]: 

/ 0 (Ara) noSSB = log (^) - 1, n = odd, (2.4) 

/(e ) noSSB | £ _ 5 . 0 -> /(e) aS y mp t = \ (log (-) +l) - (2.5) 

and in the Higgs model (1.1) Ref. [8] respectively: 

f 0 (\n) SSB = log (At)-1, (2.6) 

/(£) SSB |„0 -2 /(£)„ympt = | ( lo g((;) + l) - Y5 E ' < 2 ' 7 ) 


These results arise from integrating the known expressions [7, 8] for the tree-level amplitudes 
near the multi-particle thresholds, 


noSSB = n ! 
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over the Lorentz-invariant phase-space, a n = r Ji /' |Al ra j 2 , in the large-n non-relativistic ap¬ 
proximation. In particular, the ubiquitous factorial growth of the large-n amplitudes in (2.8)- 
(2.9) translates into the ^\A. n \ 2 ~ n!A n ~ e «iog(An) f a( q or the cross-section, which determines 
the function fo(Xn) in (2.4) and (2.6). The energy-dependence of the cross-section is dictated by 
f(s) in Eq. (2.1), and this function arises from integrating the e-dependent factors in (2.8)-(2.9) 
over the phase-space, giving rise to the small-e asymptotics in (2.5),(2.7). 

1 Equation (2.1) can be equivalently written in the form [7] a n ~ exp [A -1 e)], using the rescaling 

J-(\n,e) = XnF(Xn,e), which points towards a semi-classical interpretation of the rate in the A —> 0 limit. 
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Figure 1: Representative Feynman diagrams for the g g —> 7 h process in the SSB theory (1.1). 

An important for our forthcoming analysis point to make, is that the separability or fac¬ 
torisation of the An- from the e-dependence on the right hand side of Eq. (2.1), is the general 
consequence of the tree-level approach, i.e. it does not require taking the non-relativistic limit 
£<1 , This is because the entire A-dependence of the full tree-level result a n oc X n is contained 
in the /o function. Hence, given that the dependence on n could enter /(e) only in the combi¬ 
nation An and that the dependence on A is already fully accounted for, 2 the function /(e) does 
not depend on n. We thus should be able to determine /(e) from the fixed-n direct calculations 
of the cross-sections, and using the large-n scaling arguments suggested by (2.1). In section 3 
we will proceed to construct /(e) from the cross-sections data at n = 7 which will be computed 
numerically using MadGraph [6]. 

It should also be kept in mind that the holy grail function in the cross-section formula 
Eq. (2.1) contains not only the tree-level contributions but also the loop contributions with an 
arbitrary number of loops, which to a large extend give the dominant contributions for Xn > 1, 
as will be explored in more detail below. 


3 Results 

With MadGraph 5 [6] we can compute quite efficiently total cross-sections for scattering pro¬ 
cesses with 2 —> 7 particles. As we are mainly interested in producing the multi-particle final 
state we can make certain simplifications with respect to the characterisation of the initial 2- 
particle state. First, we take that the scattering process proceeds originates from the gluon 
fusion in the initial state, producing the highly-virtual single higgs boson h* via the effective 
ggh vertex htvG^G^, which is followed by the 1* —» n process computed in the scalar- 
theory (1.1) or (1.2), 

gg^h*^-nxh. (3-1) 

In principle, as is well-known, the use of the point-like effective ggh coupling approximation is 
not justified at high energies for producing realistic cross-sections. Finite top mass effects in loop 
will result in the form-factor in front of the exponential factor in the cross-section. However, in 
our case the single effective vertex is only a gimmick - in practice we will be computing ratios 
of the cross-sections at the same values of energy for different n. These ratios are insensitive 

2 We recall that this argument applies at tree level. 
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Figure 2: Selection of data points for log <77 as the function of e computed with MadGraph in the 
models with and without spontaneous symmetry breaking (the upper and the lower contours). 


to the bad high-energy behaviour of the effective vertex in the initial state. In our final results 
plotted in Fig. 6 we will include the effect of the Higgs form-factor, as shown in Eq. (3.12). 

At n = 7 in the Higgs theory with SSB (1.1) MadGraph computes 34,300 diagrams (two of 
which are shown in Fig. 1) contributing to the tree-level scattering amplitude. The cross-section 
values er 3SB were then computed for different energies on a grid of ~ 30 points with values of 
e = E/{7A1) — 1 ranging from 0.001 (nearly at the mult-iparticle threshold) to 250 (ultra- 
relativistic final state). We have chosen M = 125 GeV and set A = l 3 . To give an example, 
at e = 1 which corresponds to 1750 GeV, the rate is cr 3SB = 8.913 X 10 _n ± 2.74 x lCV 13 pb, 
and at e = 30 corresponding to 27125 GeV, the rate is cr 3SB = 2.818 x 10 -10 ± 9.02 x 1CT 13 pb. 
More data points for cr 3SB for 0.1 < e < 250 are shown in Fig. 2 as the upper contour (in blue). 

3.1 Extracting /(e) from computed cross-sections 

The expression in (2.3) contains only the contributions growing with n; it does not include 
subleading corrections, and is valid, as it stands, only in the n —> 00 limit. To be able to 
work at moderately large values of n, such as n = 6, 7, we now generalise this by including the 
sub-leading corrections in n. In general, they are of the form O(logn) and O(n 0 ), so that in 
total we have, 

logo/ ree = n (fo(Xn) + /(e)) + c 0 logn + ci/i(e) + c 2 + (D(l/n). (3.2) 

Here Co, c\ and c 2 are some unknown constants, and /i(e) is a new function of the kinetic energy. 
For example, by carrying out the phase-space integration of the non-relativistic amplitudes in 
the small-e limit beyond the leading order in n one finds, 

log cr/ ee ->• n (/o(An) + (3/2) log e) + c 0 logn - (5/2) log e + c 2 + 0(l/n). 

3 This can - and for the applications in plots in Figs. 5, 6 will be - rescaled to the physical value A ~ 1/8 
using the fact that for each n, o/ ree oc A” -1 . 
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Figure 3: Plots of /(e) extracted from the log CT7 ree /<7g ree MadGraph data in the SSB model, 
and the log a^ ee /a^ ee realisation of /(e) in the model without SSB. The results perfectly match 
/(s)asympt for £ < 1 depicted in light blue. 


Returning to Eq. (3.2) we now consider the difference between the rates at n and n — 1. This 
allows us to extract /(e) directly from the logcr n — logcr n _i data as follows: 


/(e) = log <Jn - fog<7n-l 


Tl 

nf 0 (Xn) - (n — l)/ 0 (A(n — 1)) + 0.5 log—— 

n — 1 


(3.3) 


The main point is that the expression in square brackets is known as it is dictated by the 
known function fo(Xn) in (2.6). (In addition, the constant cq = 0.5 is fitted from the data, and 
it results in a small correction numerically.) Equation (3.3) is our main tool for computing the 
holy grail function in the model with SSB from the ratios of cross-section data for n = 7 and 
n = 6 particles in the final state. 

The scattering amplitude into the final state with n = 6 bosons in the model (1.1) contains 
2,485 Feynman diagrams at tree level. This is still a large enough number of diagrams (to 
be in the regime of a ‘high-order’ perturbation theory), so we can use the improved large-n 
subtraction formula (3.3). The cross-sections <jg SB are computed on the same e grid as before. 
The characteristic value at e = 1 is now <jg SB = 1.77 x 10 _9 pb, and at e = 30 the rate is 
1.649 x 10~ 9 pb. 

Our results for the function /(e) in the Higgs theory (1.1) derived form the numerical cross- 
section data using (3.3) with n = 7, are shown on the left plot in Fig. 3. This plot also shows 
a perfect match to the known /(e) asympt expression (2.7) at e < 1, which is shown as a dashed 
curve in light blue. As another test of self-consistency of our procedure, we have checked that 
/(e) obtained from the 7 — 6 computation in fact matches closely the function extracted from 
a similar 6 — 5 computation. 

For completeness, and to compare with the numerical predictions based on the semi-classical 
analysis in [5, 9, 10], we have also computed /(e) in the unbroken theory (1.2). Our results 
for for cr" 0SSB for 0.1 < e < 250 are shown in Fig. 2 as the lower contour (in purple). The 
two values of the cross-sections in the two models appear to converge in the UV. This is not 
surprising, since at very high energies, all mass parameters become irrelevant and there is little 
difference between the models with the ‘right’ and the ‘wrong’ sign of the mass-squared term. 
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Figure 4: Plots of /(e) in the broken and unbroken theory for medium toi large values of e. In 
the UV regime the functions asymptote to /(e = 250) SSB ~ —2.2 and /(e = 250) noSSB ~ —1.4. 
The plot on the right depicts /(e) noSSB and shows that it exceeds the asymptotic lower limit 
— log7r 2 /2 ~ —1.6 obtained from the 0(4) symmetric classical solution [9, 10, 5]. 




To determine /(e) in the unbroken theory from the diagrammatic computation, we use the 
master formula, 


2 /(e) = log o n - logCJ n _2 


Tl 

nfo(Xn) - (n - 2)/ 0 (A(n - 2)) + 0.5 log-- 

n — 2 


(3.4) 


for the n = 7 and n = 5 rates (in the unbroken theory there are no 3-point vertices and the 
amplitudes are non-vanishing only for odd values of n). Our results are shown on the right plot 
in Fig. 3. The left panel in Fig. 4 plots the results for /(e) functions in the SSB model and the 
unbroken theory side by side for moderate to large values of e. 

In principle we should keep in mind that our analysis is based on the applicability of the 
subtraction formulae which assume that n is large enough to ensure that 1/n corrections are 
negligible in the subtraction formulae. Thus our derivation of /(e) in the unbroken theory, 
which is based on the n = 7 and n = 5 data with (280 and 10 Feynman diagrams) is less robust 
in comparison to our main SSB theory results based on the n = 7 and n = 6 data with 34,330 
and 2,485 Feynman diagrams. However computations of 2 —> 9 processes with MadGraph, 
which would be the next step in the unbroken theory, is beyond the scope of this paper. 


3.2 Multi-particle cross-sections 

Having determined the n-independent kinetic energy function /(e) allows to us to compute 
multi-particle cross-sections at any n in the large-n limit. The tree-level multi-particle cross- 
sections <7* ree in the scalar theory with SSB are obtained via 

log a n (E) = n (/o(An) + /(e)) , (3.5) 

with e(E,n) = {E — nM)/{nM ) and 

/o(An) SSBtree = log (^) 
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(3.6) 
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Figure 5: Plots of multi-particle tree-level cross-sections cr* ree in the scalar model with SSB as 
the function of energy E for a range of final-state multiplicities between n = 1000 and n = 1500. 


and we set A = 1/8 and M = 125 GeV. In Figure 5 we plot the cross-sections o// ee in this 
theory as a function of energy E for a range of final-state multiplicities between n = 1000 and 
n = 1500. The choice of such high values of particles in the final state follows from selecting the 
regime where the tree-level cross-sections become unsuppressed. This occurs when the positive 
/o(An) SSBtree factor is able to compensate the negative values of f(e). As the result we see 
that perturbative cross-sections grow very steeply with energy, and the interesting range of 
energies where the log <r* ree crosses zero occurs is the E ~ 500 TeV regime. At these energies 
the tree-level cross-sections grow exponentially violating perturbative unitarity. The energy 
regime where this happens in Fig. 5 is in agreement with the estimates obtained in [11], What 
is interesting, is that the energy scales where perturbation theory breaks down (judging from 
the leading tree-level analysis here) occurs at energies only a (few)xlO 1 above what could be 
directly tested experimentally with a hadron FCC collider. 

Let us now consider the effect of loop corrections. The 1-loop corrected multi-particle 
amplitudes on multi-particle thresholds are known [12, 13], and the result in the broken scalar 
theory (1.1) is given by [13], 

SSB : 4^ ll0 ° P = (2v) 1 ~ n ^1 + n(n - 1)^^ . (3.7) 

It was shown in Ref. [7], based on the analysis of leading singularities of the multi-loop ex¬ 
pansion around singular generating functions in scalar field theory, that the 1-loop correction 
exponentiates, 

A 1 ] 0 *^ = x exp [B Xn 2 + 0(\n)\ (3.8) 

in the limit A — > 0, n —> oo with A n 2 fixed, where B is the constant factor determined from the 
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Figure 6: Results for multi-particle cross-sections <7n'° p with the leading-loop-resummation fac¬ 
tor (3.13)-(3.14) and the single Higgs production form-facror in the model with SSB (1.1). The 
logarithm of the cross-section (3.12) is plotted as the function of energy for a range of final-state 
multiplicities between n = 110 and n = 150. 


1-loop calculation, 


model with SSB (1.1) : 

B = 

8tt ’ 


(3.9) 

unbroken model (1.2) : 

647T 2 

(log(7 + 4y/3) - m j . 

(3.10) 


As the result, the the leading-order multi-loop exponentiation leads to the the exponential 
enhancement of the multi-particle cross-section in the Higgs model, cf. Eqs. (2.6), (3.14), 

/ 0 (An) loop = log - 1 + 2B An. (3.11) 


Finally we can also include the single Higgs production form-factor in front of the exponential 
factor in the cross-section, to correct for our use of the effective Higgs-gluon vertex in the large 
energy limit, 4 



a l ° op = (mt/E) 4 \og A (rn.t/E) 2 e n (/o(An) + /(e)), 

(3.12) 

where mt is the top 

mass. In total we have 


log^°°P 

= n ^/ 0 (Ati) Ioop + /(£)) - 4 (log (E/m t ) — log log(S 2 /m,)) , 

(3.13) 

/o(An) loop 


(3.14) 


4 I would like to thank Michael Spira for this suggestion. 
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where the last equation is consistent with (3.9) and leads to the exponential enhancement of 
the cross-section o n , at least in the leading order in n 2 A. The form-factor correction - the 
last term on the right hand side of (3.13) - grows with n only logarithmically 5 compared to 
the linear in n terms in the first term. (At E = 50 TeV the form-factor gives the correction 
— 4 (log (E/mt) — log log(E 2 /m 2 )) ~ —12.95 in the exponent.) 

Our results for a^n ° P including the form-factor and the exponentiated loop factor (the last 
term in (3.14)) for the Higgs model (1.1) are shown in Fig. 6 for a range of final-state multi¬ 
plicities between n = 110 and n = 150. We can see that the loop-enhancement has reduced 
the energy scale (and multiplicities) by a factor of 10, and the scale for the perturbation theory 
breakdown or equivalently the scale of new phenomena to set in is now reduced to 40-50 TeV 
which is pretty much within the energy reach of the 100 TeV FCC collider, in agreement with 
the estimate in [11]. 

Of course, one should keep in mind that the set-up in Eqs. (3.13)-(3.14) is merely an opti¬ 
mistic phenomenological model. In general the even higher-order effects of loop exponentiation 
will be present such that, 

/o(An) allloops = log - 1 + ^3 ^ + const + const' + ... , (3.15) 

and can change the cross-sections contours in Fig. 6. (Note that the value of the loop expansion 
parameter are ~ 1 for n = 100 and ~ 1.4 for n = 140.) 


4 Conclusions 


Our diagrammatic approach is conceptually different (but also complimentary) to the semi- 
classical considerations followed in the earlier literature. The exponential form of the cross- 
section in the large-n limit, 


o n ~ exp 


1 

A 


AnF(An, e) 


exp 


1 

A 


F(\n, e) , 


(4.1) 


is strongly suggestive of an underlying semi-classical origin of the multi-particle cross-section. In 
particular, there is a strong similarity between the purely perturbative multi-particle processes 
considered here and the B + L-violating non-perturbative reactions in the instanton sector of 
the Standard Model discussed originally in [14], 

The idea that semi-classical methods can be also used in the perturbative sector of the 
theory was put forward and explored by a number of authors including Refs. [9, 10, 5, 15, 16]. 
At tree level the holy grail function F or T can indeed be reconstructed numerically if one can 
determine certain singular classical solutions to the boundary value problem [9], as explained in 
[5, 9, 10]. In practice this procedure was carried out in the case of the unbroken /> 4 theory (1.2) 
and based on finding numerically the singular solutions with the hypothesised 0(4) symmetry. 
In this approach a lower bound on the tree-level cross-section (2.1) was derived in [10, 5] which 
corresponds to an upper bound on the absolute value of |/(e)|. In particular it was found that at 
infinite energies, e —> oo the function /(e) —> — log(7r 2 /2) ~ —1.6. In our case, the asymptotic 

5 At large n and large e limit it is ~ —41og(ne) + 41oglog(ne) + 4 log (( mt/M ) log (M/rrit)) ■ 
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value appears to be smaller in magnitude in the non-SSB theory, /(e) —> —1.6, as can be seen 
from the right panel in Fig. 4. The fact that /(e) —» — |const| implies, 

^treenoSSB > g - |const| n e n/ 0 (An) a t E -> OO . (4.2) 

An often quoted misreading of this result is the statement that perturbative cross-sections 
remain unobservable in the multi-particle limit, even at infinitely high energies, due to a rising 
with n exponential suppression factor e -l const l n _ This, of course, is not the case as the plots 
in Figs. 5 and 6 demonstrate: the growing function e n 7o(An) compensates the suppression in 
e~ l const l n for any e already at moderately high values of An. 

One advantage of the diagrammatic approach followed in this paper is its simplicity, and also 
the fact that one should be able to apply it in any theory, ultimately including the full Gauge- 
Higgs theory of the Standard Model weak sector by generalising the non-relativistic results of 
[17, 8] to the general-c case. We leave this to future work. 

Since in our case the calculations are carried out within the first principles perturbative 
approach, we also know that as soon as the regime is reached where the theory breaks down 
and violates unitarity, this implies that we really are falsifying the perturbative technique itself, 
and not a bound arising from a semi-classical treatment. The perturbation theory break-down 
found here occurs in two cases: a) within the tree-level approximation in the energy-multiplicity 
regime of Fig. 5, and b) within the leading-order in the loop expansion approximation in the 
regime corresponding to Fig. 6. 

The main technical challenge which still needs to be addressed is how to account for all 
the remaining higher-loop corrections, relevant in the regime An ~ 1. Even the leading-order 
exponentiation of the loop corrections result (the last term on the right hand side of (3.13)) 
which was essential for lowering the characteristic energy scale from Fig. 5 to Fig. 6 by an order 
of magnitude, has been derived only in the multiparticle threshold limit; the fool e-dependence 
of leading loop corrections remains unknown. 

Another (perhaps less crucial) technical limitation of our simple derivation is that we have 
concentrated only on sub-processes with a single virtual Higgs in the s-channel in (3.1). We 
have not considered here the effect of possible numerical partial cancellations between the s- 
channel t —> h* —> hh and box diagram processes t —> hh and in the double Higgs case and 
generalisations for the multi-Higgs case, see e.g. the discussion in [18]. However, in the absence 
of the symmetry reason, we do not expect that such partial cancellations could significantly 
modify the exponential growth of the s-channel processes (3.1). This conclusion is also n 
agreement with the discussion in section 4 of [11] where we have seen that the exponential 
growth persists in the similar case of the weak vector boson fusion, VV —>■ h* —> hh —> n x h 
vs VV —> h*h* —>• n x h. 

We have shown that in very high energy scattering events, perturbative rates for production 
of multiple Higgs bosons grow with increasing energy, eventually violating perturbative unitarity 
and resulting in the breakdown of the ordinary weakly-coupled perturbation theory. The energy 
scales where electro-weak processes can enter this regime are potentially within the reach of the 
100 TeV future hadron colliders, or at least not much above it. It was argued in [11] that novel 
physics phenomena must set in before these energies are reached: either the electroweak sector 
becomes non-perturbative in this regime, or additional physics beyond the Standard Model 
might be needed. 
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